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Abstract 

The scalar radius of the pion plays an important role in xPT, because it 
is related to one of the basic effective coupling constants, viz. the one which 
controls the quark mass dependence of at one loop. In a recent paper, 
Yndurain derives a robust lower hound for this radius, which disagrees with 
earlier determinations. We show that such a bound does not exist: the 
"derivation" relies on an incorrect claim. Moreover, we discuss the physics 
of the form factors associated with the operators uu, dd and ss and show 
that their structure in the vicinity of the KK threshold is quite different. 
Finally, we draw attention to the fact that the new data on the slope of 
the scalar K^^ form factor confirm a recent, remarkably sharp theoretical 
prediction. 



1. Introduction 



Early work on the scalar form factors of the pion jTjE] was motivated by the search 
for a very light Higgs particle. Unfortunately, the outcome of this search was 
negative: nature is kind enough to let us probe the vector and axial currents, but 
allows us to experimentally explore only those scalar and pseudoscalar currents 
that are connected with flavour symmetry breaking. In particular, there is no 
handle on the matrix element^ 

T^(t) = {7i{p')\muuu + md(Id\7i{p)) , t = {p' - pf . (1) 

The value of this form factor at t = is referred to as the pion cx-term, 

dMl dMl 
orriu orrid 

According to Gell-Mann, Oakes and Renner |^, the expansion of the square of 
the pion mass starts with a term linear in and m^, the coefficient being 
determined by the quark condensate. Hence the pion a-term is determined by 
the pion mass, except for corrections of higher order: r7r(0) = + O(m^), with 
rh = + md)- Indeed, both the precision measurements on Ki4^ decay and 
the preliminary results from DIRAC 5J confirm that the corrections are small: 
more than 94 % of the pion mass originates in the term generated by the quark 
condensate p. 

The scalar radius represents the slope of the corresponding normalized form 

factor f^(t) = r^(t)/r^(o), 

TAt) = l+^-{r'):t + Oit'), (3) 

which is of considerable interest, because it is related to the effective coupling con- 
stant £4, that determines the first nonleading contribution in the chiral expansion 
of the pion decay constant. Denoting the value of in the limit = = 
by F, we have [2] 

There is a formula analogous to Q also for Fk/ F^^. Neglecting Zweig rule vio- 
lating contributions and using the measured value of Fk/F.,^, this relation leads 
to a first crude estimate for the scalar radius: = 0.55 ± 0.15 fm^ jSj. An im- 

proved estimate was obtained long ago on the basis of dispersion theory j2] . The 
calculation relied on the assumption that only the transition tttt KK generates 
inelasticity at low energies - all other inelastic channels in the Mushkhelishvili- 
Omnes (M-0) representation of the form factor were neglected. Moussallam 

"'^We work in the limit iriu = md, e = 0, where isospin is an exact symmetry. 
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[H] performed a thorough analysis of this approach, considering several different 
phase shift representations (in particular also the parametrizations proposed in 
|10j ) and studying the sensitivity of the outcome to other inelastic reactions, such 
as TTTT — An, TTTT — > r]7]. His results for the scalar radius are in the range from 0.58 
to 0.65 fm^. In [TT], the Roy equations for tttt scattering were used to update the 
calculation described in [2], with the result 

(r^)- = 0.61 ±0.04 fm^ (5) 

The central value confirms the number r^(0)/r7r(0) = 2.6 GeV^^ given in |2| and 
the error bar covers the range found by Moussallam. 

The higher orders of the chiral perturbation series for the form factor T^^ (t) are 
discussed in [T2] and a detailed comparison with the dispersive representation can 
also be found there. The complete evaluation to two loops is given in ^^l- The 
corrections of 0{rh?) in the relation (jH) are discussed in ^Mj. With the estimates 
for the higher order terms given in ^3], we obtain F^^/F = 1.072 ± 0.004. As the 
corrections are (a) very small and (b) dominated by known double logarithms, 
the uncertainty in the result is due almost exclusively to the one in the scalar 
radius. 

The effective couplings relevant for the masses and decay constants can be 
measured on the lattice (TH]. Using the values for L4, L5, Lq and Lg found by the 
MILC collaboration [T7j, the corresponding values of the SU(2) xSU(2) coupling 
constants are readily worked out from the relations given in P!H]. This leads to 
£3 = 0.8 ± 2.3, £4 = 4.0 ± 0.6, in good agreement with the estimates given 20 
years ago. Inserting this value for £4 in the relevant one loop formulae, we obtain 
F„/F = 1.06 ± 0.01, (r^)^ = 0.5 ± 0.1 fm^ A direct determination of the ratio 
Ft,/F on the lattice would be of considerable interest. 



2. Omnes representation 

Yndurain's paper on the subject [1^] is based on the Omnes representation, 

t =exp- / , 6 

TT 74M2 S{S — t) 

which expresses the form factor in terms of its phase on the upper rim of the cut, 
6y{s) = argr7r(s + ie). The formula may be viewed as a one-channel version of 
the M-0 representation (in that framework, the absence of inelastic channels im- 
plies that the phase of the form factor coincides with the phase of the scattering 
amplitude). Perturbative QCD indicates that the form factor behaves asymptot- 
ically as |r^(t)| ~ l/|t| up to logarithms [201 • form factor does not have 
zeros, the phase Sr{s) must tend to vr. The formula © then rigorously holds and 
leads to a rapidly convergent representation for the scalar radius, 

6.»<w^^ (7) 

^ TT J AMI ^ ' 



3 



Unless the asymptotics is assumed to set in at an unreasonably low energy, the 
corrections from the preasymptotic logarithms are negligibly small. 

The Watson theorem states that, in the elastic region, the phase of the form 
factor coincides with the isoscalar S-wave tttt phase shift, 6r = S^. Below the 
KK threshold, inelastic processes do not play a significant role: on the interval 
4M^ < s < 4M|^, the elasticity 7]q remains very close to 1, so that 6r remains 
close to 6q. With the representation for the phase 5q obtained in ref. pQ], the 
contribution from that interval can be evaluated quite accurately. 

The opening of the KK channel produces a square root singularity at 4M|-, 
which manifests itself as a dip in the elasticity in the region between 1 and 1.1 
GeV. Although the valley may not be very deep, there is no reason for the phase 
of the form factor to agree with 6q in that region. Various models have been 
proposed to account for the fact that the Omnes factor belonging to Sq does not 
properly describe the behaviour of the form factor r^(t) or of other transition 
amplitudes involving the production of pion pairs (see for instance PHI^I^ I24j). 

Yndurain assumes that the perturbative asymptotics sets in at 1.42 GeV, ob- 
serves that in the region between 1.1 and 1.42 GeV, the inelasticity is compatible 
with zero and then claims: It thus follows that the phase o/r^(s) must be approxi- 
mately equal to 5q{s) for 1.1 GeV < s^/^ < 1.42 GeV. This claim is incorrect, for 
the following reason: in the presence of inelastic channels, the Watson theorem 
in general reads 

T*^is)=Y.{5mn + 2lT^nis)anis)yTr,{s). (8) 

n 

We use the notation of |2I and identify the first two channels with vrvr and KK: 

2 

ri(s) = r^(s), r2{s) = --^Tk{s) , o-i(s) = o-^(s), 0-2(3) = o-i^(s), (9) 

with ap{s) = 6{s - 4M|,)(1 - 4M|,/s)^/^ The term Tn stands for the partial 
wave amplitude of the isoscalar S-wave, 

r/gexp(2z^g)-l 
= ^° 2^^^ • 

If all other channels are ignored, unitarity fixes the magnitude of T12 above the 
KK threshold in terms of the elasticity: 4(Ticr2 |Ti2p = 1 — iVo)"^- For energies 
where rj^ ~ 1, the condition (jH)) thus reduces to F* ~ exp(— 2 iS^) F^. This relation 
does not imply that the difference 6^ — approximately vanishes, but only 
requires that it is close to a multiple of vr. 

One might think that continuity would remove the ambiguity, but this is not 
the case, because the region of interest is separated from the elastic domain by 
an interval were inelasticity cannot be ignored. The full line in fig. ^ depicts the 
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Figure 1: The full line is the phase of the pion form factor of the operators 
uu or dd, as calculated from the two-channel M-0 equations. The dashed and 
dotted lines describe the corresponding phase shift and the phase of the partial 
wave amplitude, respectively. The dash-dotted line depicts the phase of the form 
factor belonging to the operator ss. 

outcome of our calculation^ for the phase of the form factor: above 1.1 GeV, 6r 
indeed differs from 6q, approximately by vr. The detailed behaviour in the region 
around 1 GeV is sensitive to the properties of the T-matrix, but the entire range 
of representations considered in ref. P leads to a sharp drop of Sr at the KK 
threshold and to (5q — (5r — vr for energies above 1.1 GeV [26] . In other words, 
the robust lower bound in [1^] is not valid, because it is derived from an incorrect 
claim. 

The dotted line in the figure depicts the phase of the partial wave amplitude, 
6t = arg^Q - the phase of the form factor closely follows this line. The explicit 
calculation based on two-channel unitarity thus leads to a behaviour of the scalar 
form factor of the type proposed by Morgan and Pennington for the diffractive 
production of vrvr final states [21]. Indeed, fig. 2 in their paper on the reaction 
pp — >• ppmi j2Zl is closely related to ours: it shows the result of the energy 

^The specific curves shown in the figure are based on the T-matrix representation of Hyams 
et al. [53] . More precisely: (a) that representation is used as it is only on the interval 0.8 GeV < 
E < 1.5 GeV; (b) at lower energies, we fix Tn as well as the phase of T12 with the solution of 
the Roy equations specified in (17.1), (17.2) of ref. JT], taking only the ratio |ri2/T'ii| from the 
Hyams representation; (c) on the interval from 1.5 to 1.7 GeV, T is guided to zero smoothly, 
in accordance with the unitarity condition {6q ^ 2 7r, argip tt). 
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independent analysis of the Hyams data for Sq and 6t, while, above 4M|', the 
dotted and dashed curves in our plot represent the result of the energy dependent 
analysis of the same data. 

3. Importance of inelastic channels 

The reason for the pronounced difference between 5° and 5t is readily understood 
from the Argand diagram: As the energy reaches 2Mk, the amplitude has nearly 
completed a full circle. If inelasticity could be ignored, the curve would continue 
following the Argand circle, so that would have a zero, a few MeV above the 
KK threshold. Hence the phase 6t would make a jump there, dropping abruptly 
by vr. In reality, the curve leaves the circle before the phase has reached vr, so 
that t[] remains different from zero and a jump does not occur. Instead, 6q — 6t 
continuously, but rapidly grows from to the vicinity of vr. 

The phenomenon illustrates the fact that phases of small quantities can be 
very sensitive to details. The phase of undergoes a dramatic change because 
it so happens that an inelastic channel opens up at an energy where nearly 
vanishes. The magnitude of the change in 5^ — 6r is by no means proportional to 
the probability for the formation of a KK pair, i.e. to the inelasticity 1 — (?7o)^5 
but is approximately equal to vr. If the inelasticity is small, the change in the 
phase difference takes place almost instantly. 

In connection with the Omnes formula, the difference between the phase shift 
and the phase of the partial wave is a measure of the importance of inelastic 
channels. Above 1.1 GeV, both — and S-p — obey the Watson theorem. 
Fig. 121 shows that in the region below 1.4 GeV, an evaluation of the integral rele- 
vant for the scalar radius based on 6r — St practically reproduces the result of our 
two-channel calculation, while using 6r — leads to values hke those advocated 
in ^^1; which are significantly higher. So, inelastic reactions are important here: 
In order to determine the scalar radius, we need to know their impact on the 
form factor. 

For the electromagnetic form factor of the pion, the situation is qualitatively 
different. In that case, inelastic channels play a much less important role. In 
particular, the angular momentum barrier suppresses the branch point singularity 
connected with the opening of the KK channel. Since the P-wave phase shift Si 
stays well below 180°, the partial wave amplitude t\ does not become small there, 
so that the phenomenon observed in the S'-wave does not occur: the difference 
between 6l and the phase of t\ grows much more slowly than in that case: at 
1.1 GeV, it amounts to a few degrees, while 6q ~ 6t — 180°. In this connection, 
we recall that, for the case of the electromagnetic form factor, Eidelman and 
Lukaszuk have shown that the experimental information on e~^e~ production 
of final states other than tttt implies rather stringent bounds on the elasticity rjl 
and on the difference between the phase shift Si and the phase of the form factor. 

The behaviour of the form factor at the onset of the KK continuum reflects the 
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strength of the couphng to these states. For the operator ss, this couphng differs 
from Tiu or dd. Hence we should expect that the phase (5a (s) = arg A,r(s + « e) of 
the form factor 

A^{t) = {7r{p')\msSs\7T{p)) (11) 

behaves quite differently from 5r('S). As shown in |2|, Ajr(t) is given by a different 
linear combination of the same two linearly independent solutions of the M-0 
equations that are needed for the evaluation of r^(i(:). In fig. [H the phase of the 
resulting representation of the form factor is shown as a dash-dotted line. The 
figure shows that roughly follows the phase shift 6^: above the KK threshold, 
the phases of the two form factors are very different. 

4. Magnitude of the form factors 

Fig. 121 shows the magnitude of the Omnes factors obtained by inserting the phases 
depicted in fig. Q in the formula (jH)). The full curve represents our result for 
the form factor |r^(s)| and shows that this quantity exhibits a dip at the KK 
threshold - the phenomenon discussed by Morgan and Pennington. Indeed, the 
figure shows that the result for this form factor is very close to the Omnes factor 
belonging to the phase 6t. Moussallam's analysis confirms the phenomenon: for 
all of the T-matrix representations considered in |9j, the function |r^(s)| goes 
through a sharp minimum near the KK threshold 

The minimum reflects the rapid drop in the phase: the Omnes factor belonging 
to the phase 6{s) = 6{s — 4M|-) (— vr) is given by 1 — t/AMj^. In other words, 
if the phase were to drop suddenly by tt at s = 4M|^, then the corresponding 
Omnes factor would contain a zero there. In reality, the phase does not drop 
suddenly, but rapidly - the magnitude of the form factor does not go through a 
zero, but through a minimum. Conversely, the fact that the form factor becomes 
very small near the KK threshold implies that the behaviour of its phase there is 
very sensitive to details and cannot be understood without explicitly accounting 
for the KK channel. 

For 5a ; on the other hand, the Omnes factor exhibits a peak near the KK 
threshold. Below 1 GeV, the behaviour is very similar to the one of the Omnes 
factor evaluated with 6^: if (as advocated by Yndurain) the phase of were to 
follow 6q rather than 6t, this form factor would exhibit a pronounced peak rather 
than a dip. This reflects the fact that the operator ss couples more strongly 
to the kaon than to the pion. Near t = 0, the difference is not enormous, but 
the slope is of course larger: evaluating the integral in (0) with (5a instead of 
Sr, we obtain 0.81 fm^, instead of the number 0.61 ± 0.04 fm^ quoted above. The 
behaviour of 6a near the KK threshold is subject to considerable uncertainties 
- we did not make an attempt at estimating those in the corresponding radius. 

There is a qualitative difference between the two form factors under consid- 
eration here: our representation for A^^it) has a zero, but F^(t) does not. The 
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Figure 2: Omnes factors belonging to the phases shown in fig. 1. 

reason is that At,{0) represents the derivative of with respect to and hence 
vanishes for m„ = rrid = 0, while the slope A^(0) does not disappear in that limit. 
Accordingly, the Omnes representation involves a polynomial: 

A. t = Po + Pit exp- / 12 

TT J4A/2 S{S — t) 

The explicit representation in xPT to one loop jH] shows that po is of 0{m), while 
Pi is of 0(1) (the representation exclusively involves the Zweig rule violating 
constants L4 and Lg). This demonstrates that the form factor A-j^lt) necessarily 
has a zero at a value of t of order rh, i.e. in the region where xPT is reliable. In 
order for the representation (fT^ to be consistent with perturbative asymptotics, 
the phase 6a must tend to 2tt (compare fig. H}. Note that the dash-dotted curve 
in fig. 121 represents the magnitude of the exponential and does not account for 
the polynomial. 

5. Scalar radius relevant for Ki^ decay 

The scalar form factor relevant for the decay K — > vrfz/ is proportional to the 
matrix element {K\ su\tc). We denote this form factor by fo{t), using the standard 
normalization, where the value at t = coincides with /+(0), a quantity that is 
of central importance for the determination of the CKM matrix element Vus- The 
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term linear in t, 

m = UiO) [l + l {r')f-t + Cot' + ...y (13) 

represents an analogue of the scalar radius of the pion. In the analysis of the data, 
it is customary to replace this radius by the slope parameter Aq = (r^)f'^M^/6. 

Nearly 20 years ago, a prediction for the radius was made, on the basis of xPT 
to one loop: (r2)f^ = 0.20 ± 0.05 fm^ |H]. This numb er is about 3 times smaller 
than the scalar radius of the pion, an illustration of the fact that the scalar radii 
are very sensitive to flavour symmetry breaking - in contrast to the vector radii, 
where the flavour asymmetries are comparatively small. The corrections to the 
Callan-Treiman relation were also analyzed. In the formulation of Dashen and 
Weinstein, this relation represents a low energy theorem ^j, which states that 
in the limit = = 0, the value of fo{t) at t = M"^ — coincides with the 
ratio Fk/Ft^. As it turns out that the corrections of 0{rh) do not contain a chiral 
logarithm of the type logM^, they are tiny P^. 

The experimental situation was not clear at that time: The outcome of a high 
statistics experiment jHO] was in agreement with the theoretical expectations, but 
as explicitly stated in [S], the values for (r^)f"'^ found in some of the more recent 
experiments cannot be reconciled with chiral symmetry. In the analysis of the 
Particle Data Group, the unsatisfactory experimental situation manifests itself in 
the fact that (a) the scale factors S needed to account for the inconsistencies are 
large and (b) despite the stretching of error bars, the value found from decays of 
neutral kaons does not agree with the one from decay. 

In this field, there was considerable progress recently, on the experimental 
as well as on the theoretical side. In particular, the Ki^ form factors are now 
known to two loops of xPT pT| l32j. The curvature of the form factors cannot be 
neglected at the precision reached now and, in principle, a precise experimental 
determination thereof would allow a parameter free measurement of Vus |32j . 
Moreover, Jamin, Oiler and Pich observed that the curvature of the scalar form 
factor can be determined rather accurately by means of dispersive methods j33j . 
This implies that the radius can be calculated from the value of /o(t) at the 
Callan-Treiman point, t = M]^ — M^, for which xPT makes a very accurate 
prediction. In this way, the authors arrive at 

(r2)f^ = 0.192 ±0.012 fm^ ISni. (14) 

The central value confirms the old result mentioned above, the uncertainty is four 
times smaller. 

In p!^, Yndurain states that the value of Aq for charged kaon decay published 
by the PDG in 2000 is difficult to believe. Discarding the data prior to 1975, he 
arrives at {r"^)^^ = 0.312 ± 0.070 fm^ and concludes that the central value lies 
clearly outside the error bars of the chiral theory prediction. Indeed, if his central 
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value was close to reality, we would have to conclude that experiment is in flat 
contradiction with a low energy theorem of SU(2) xSU(2). 

This is not the case, however. For the charged kaons, the data collected at 
the ISTRA detector clarified the situation considerably ^Mi- The result for the 
radius reads (r^)f = 0.235 ± 0.014 ± 0.007 fm^ (note that in this case the radius 
is calculated from Aq using M^^o) , which now dominates the world average. For 
the neutral kaons, the experimental situation also improved significantly: there 
IS a new result from KTeV, (r^)^^^ = 0.165 ± 0.016 fm^ Since this value (a) 
now dominates the statistics and (b) is consistent with the 1974 high statistics 
experiment mentioned above, we conclude that there is a problem with those 
of the earlier data that were in conflict with chiral symmetry. While the value 
obtained from decay is higher than the prediction ()14j) by 2.1 cr, the KTeV 
result is lower by 1.4 a. Chiral symmetry indicates that the truth is in the middle. 

8. Conclusion 

1. The low energy properties of the scalar pion form factors are governed by those 
of the isoscalar S-wave in tth scattering. In particular, the reaction tttt — >■ KK 
generates a pronounced structure in the vicinity of s = 4M|., which can be 
understood on the basis of a dispersive two-channel analysis. This framework 
leads to the conclusion that, in the region around 1 GeV, the pion matrix elements 
of uu and (Id roughly follow the vrvr partial wave amplitude and thus exhibit a 
sharp minimum there. The coupling of the operator ss to the KK states differs 
from the one of uu or dd. The corresponding form factor exhibits a pronounced 
peak rather than a dip. 

2. The dispersive analysis leads to a rather accurate determination of the 
scalar radius of the pion. The early estimate given in |2| is confirmed. In par- 
ticular, as shown in ^1, the uncertainties in the phenomenological information 
used above 0.8 GeV do not significantly affect the result, which is in the range 
(r^)g = 0.61 ± 0.04 fm^ jTTj. We draw attention to the fact that the two-loop 
prediction for the dependence of the pion decay constant on the mass of the two 
lightest quarks jHj can be used to convert determinations of on the lattice 
into a measurement of the scalar radius. The existing lattice data are consistent 
with the result of the dispersive calculation. 

3. We discuss the impact of the new precision data on the scalar form factor 
of Kg3 decay . Chiral symmetry leads to a low energy theorem for the 
value of this form factor at t = M"^ — M^. The new results, which now dominate 
the statistics, show that there is a problem with those of the old data that were 
in conflict with this prediction. Combining xPT to two loops pTl with a 
dispersive analysis of the curvature, the low energy theorem can be converted 
into a very sharp prediction for the radius {r'^)f'^ or for the slope parameter Aq 
PH] . Unfortunately, in view of the very small errors quoted for the slope, the new 
data on Ki decay are not in agreement with those on K^ decay: while the former 
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are lower than the prediction, the latter are higher. Hopefully, the analysis of 
the data collected by the KLOE collaboration at Frascati and by NA48 at 
CERN inZj will clarify the situation. 

4. Yndurain jlH] states that the two-channel analysis in ^ is of the "black- 
box" type and claims that it is not necessary, that the phase of the form factor 
must approximately follow the phase shift Sq, that the scalar radius of the pion is 
subject to a lower bound and that the chiral theory prediction for (r^)f disagrees 
with experiment. We have shown that none of these claims is tenable. 
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